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Recent results from lattice QCD calculations relevant to particle physics phenomenology are reviewed. 
They include the calculations of strong coupling constant, quark masses, kaon matrix elements, and D 
and B meson matrix elements. Special emphasis is on the recent progress in the simulations including 
dynamical quarks. 



1 Introduction 

Since it was invented by K. Wilson 1 in 1974, 
lattice QCD has grown into an important 
tool for the analysis of low energy regime of 
QCD, where the non-perturbative dynamics 
of quarks and gluons becomes essential. In 
this talk I review the status of the calcula- 
tion of several quantities relevant to particle 
physics phenomenology, such as the determi- 
nation of fundamental parameters of QCD, 
and the hadron matrix elements of K, D and 
B mesons. 

Although lattice gauge theory starts with 
first principles in QCD calculations, due to 
the limitation of computational resources one 
has to resort to several approximations, such 
as finite lattice spacing a, finite spatial extent 
L of the lattice, and relatively large quark 
mass m q . They introduce sources of system- 
atic uncertainties, which can in principle be 
reduced by performing extrapolations to the 
appropriate limits. 

Among other approximations, the most 
serious one was the quenched approximation, 
in which gluon vacuum polarization of by 
quarks is neglected. It reduces the compu- 
tational cost by orders of magnitude and so 
has been used for most lattice calculations, 
but the systematic uncertainty is not un- 
der any quantitative control. The most im- 
portant development since the last ICHEP 
conference in 2002, where Lellouch summa- 
rized the lattice results 2 , is that the simula- 
tions beyond the quenched approximation — 



unqucnched QCD — have been carried out for 
many important quantities. 

In this review, I firstdiscuss the theoret- 
ical issues in dynamical quark simulations, 
which are related to formulations and algo- 
rithms for fermions on the lattice. Then I 
cover several applications, including the de- 
termination of the strong coupling and quark 
masses, and calculations of weak matrix ele- 
ments. Other reviews of recent results can be 
found in 3,4 . 

2 Issues in dynamical QCD 

Since dynamical fermion simulations involve 
many inversions of fermion matrix, it is sub- 
stantially harder to simulate light quarks at 
their physical mass values. Therefore, an ex- 
trapolation, called the chiral extrapolation, 
in the light quark mass from feasible quark 
masses to the physical up and down quark 
masses is necessary. The extrapolation is best 
controlled with the functional form predicted 
by the chiral perturbation theory (ChPT), 
which provides a systematic expansion of low 
energy QCD for small quark masses. 

The question is then whether the re- 
gion of quark masses in lattice simulations 
has enough overlap with the convergence ra- 
dius of ChPT. For the pion decay constant, 
for instance, at NLO ChPT predicts a non- 
analytic behavior — m\ mm 2 , called the chi- 
ral log. The result from the JLQCD col- 
laboration 5,6 clearly shows that there is 
no signal of curvature for the "pion" mass 
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above 550 MeV, which is currently the low- 
est available pion mass with (improved) Wil- 
son fcrmion formulations. More recently 
the MILC collaboration 7 has performed a 
simulation with the pion mass as low as 
250 MeV using the (improved) staggered 
fermions. They found that the chiral log 
shows up below 500 MeV, showing a clear ad- 
vantage of the staggered fermion. In fact, us- 
ing the simulations including 2+1 (up, down 
and strange) flavors of staggered quarks, it 
has been demonstrated that several funda- 
mental physical quantities are in good agree- 
ment with the experiment 8 . 

Wilson fermions explicitly violate chiral 
symmetry at finite lattice spacing, and the 
massless quark limit has to be reached by 
tuning a mass parameter. Because the mass- 
less point itself fluctuates statistically in the 
Monte Carlo simulation, the singularity of 
the massless limit could show up earlier in 
the chiral extrapolation, and the dynami- 
cal fermion simulation becomes much harder 
near the chiral limit than expected from naive 
scaling law. With the staggered fermion, on 
the other hand, there is an exact U(l) chiral 
symmetry and the massless limit is fixed. 

The price one has to pay for the stag- 
gered fermion is the complication of species 
doubling. The formulation necessarily in- 
volves four species (or tastes) of fermions. 
They mix among themselves at finite lattice 
spacing. As a result, there are 16 pions on the 
lattice, only one of which becomes massless 
in the chiral limit. All other hadrons are also 
duplicated, and such effects could become a 
source of systematic uncertainty. 

More seriously, one has to take a square- 
root or fourth-root of the fermion determi- 
nant in order to represent degenerate up and 
down quarks or slightly heavier strange quark 
using the four-taste staggered fermion. Any 
local Dirac operator corresponding to the 
fourth-rooted staggered fermion determinant 
has not been found so far a . Without the 



locality one cannot prove the universality, 
i.e. the continuum limit is the QCD. There- 
fore, it is potentially a fundamental problem 
of the calculations employing the staggered 
fermions 6 . 

Theoretically, the best solution is to cm- 
ploy the fermions satisfying the Ginsparg- 
Wilson relation, which have an exact chi- 
ral symmetry at finite lattice spacing with- 
out sacrificing the flavor symmetry 13 . Neu- 
berger (overlap) fermions 14,15 and domain- 
wall fermions 16 ' 17,18 fall into this class. 
Exploratory dynamical simulations using 
the domain-wall fermion have already been 
performed 19 . 

In the following discussion I assume that 
the fourth-rooted staggered fermion is a cor- 
rect, or at least effective, description of QCD 
and do not quote errors associated with it. 

3 Fundamental QCD parameters 

3.1 Strong coupling constant 

The strong coupling constant a s (Mz) can 
be determined using lattice QCD by con- 
verting o4°*(l/a) to the continuum defini- 
tion af s (Mz) by perturbation theory. The 
quarkonium spectrum is often used to set 
1/a, because it is insensitive to other sys- 
tematic errors, such as those from finite vol- 
ume and chiral extrapolation of light quarks. 
To improve the perturbative expansion the 
renormalized coupling 20 ay (<?*), which is de- 
fined through the heavy quark potential, is 
used with an appropriate scale q* . 

Figure 1 summarizes lattice results us- 
ing unquenched simulations. After includ- 
ing the dynamical quark effects, the most 
important source of systematic error is the 
unknown higher order perturbation theory. 
Previous results used two-loop matching, but 
this year the HPQCD collaboration 26 has 



a It was shown that the square-root of the staggered 



operator is non-local 9 ' 10 . 

'There is also a positive indication, i.e. the eigen- 
value spectrum of the staggered operator shows an 
approximate four-fold degeneracy 11 ' 12 . 
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Figure 1. Strong coupling constant a s (Mz)- Lat- 
tice data are from Aoki et al. 21 , Davies et al. 22 , 
SESAM 23 , QCDSF-UKQCD 24 and its update 25 , 
HPQCD-UKQCD-MILC-Fermilab 8 and its 
update 26 . 



carried out three-loop calculations and re- 
ported a (preliminary) result with a substan- 
tially reduced error, a s (Mz) = 0.1175(15), 
which is in good agreement with the PDG 
2004 average 0.1182(20) 27 . They use the sim- 
ulation data with 2+1 flavors of improved 
staggered fermions at three lattice spacings 
and confirm that various input quantities 
to determine ay(q*) give a consistent result 
within estimated four-loop errors. It should 
be noted that the results with the Wilson- 
type fermion 24 ' 25 is significantly lower. It is 
likely an unknown higher order effect, but is 
not fully understood. 

3.2 Light quark masses 

The light (up, down and strange) quark 
masses are determined from the pion and 
kaon masses, e.g., using the PCAC relation 
such as m 2 K = B(fh + m s ) at the leading or- 
der of m q . rh denotes an average up and down 
quark mass and m s is the strange quark mass. 
Since the effect of chiral log is not significant 
for m s , most of the calculations use a linear 
fit in an average light quark mass for interpo- 
lation. For the calculation of m, on the other 
hand, the NLO effect of ChPT is important. 

Because m q is regularization dependent, 
the lattice results are usually quoted in the 
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Figure 2. Strange quark mass m s (2 GeV) (MeV) 
from lattice QCD. Both quenched and unquenched 
results are listed. Quenched results (upper panel) are 
from Gimenez et al. 30 , Becirevic et al. 31 , JLQCD 32 , 
ALPHA-UKQCD 33 , QCDSF 34 , Becirevic et al. 35 , 
CP-PACS 36 ' 37 , SPQcdR 38 , Hernandez et al. 39 , 
Giusti et al. 40 , Chiu-Hsieh 41 ' 42 , DeGrand 43 , Blum 
et al. 44 , CP-PACS 45 , RBC 46 . Unquenched results 
are from CP-PACS 36 , JLQCD 6 , QCDSF-UKQCD 47 , 
SPQcdR 48 , CP-PACS/JLQCD 49 , HPQCD-MILC- 
UKQCD 50 - 7 . PDG 2004 27 average is shown by a 
dashed band. 



MS scheme by using perturbative matching 
and sometimes also using non-perturbative 
techniques at intermediate steps. For the 
case of the Wilson- type fermions, the deter- 
mination through axial Ward identity (AWI) 
and vector Ward identity (VWI) could be dif- 
ferent at finite lattice spacing c . 

In Figure 2, I compile the lattice results 
for strange quark mass for both quenched and 
unquenched calculations. In the quenched 
approximation, systematic studies of the non- 
perturbative matching and the continuum 
extrapolation have been extensively studied 
and the results are in agreement within sa 
10% quenching error, which appears as a de- 
pendence on the input quantity to set the lat- 
tice scale, e.g. m p , fx, etc. 

The CP-PACS 36 and JLQCD 6 collabo- 



c For further discussions on the quark mass calcula- 
tions I refer to 28 - 29 . 
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rations found that m s (2 GeV) becomes sig- 
nificantly lower by the effects of two dynami- 
cal flavors. d Recently, 2+1-flavor calculations 
are reported by the CP-PACS/JLQCD 49 
and HPQCD-MILC-UKQCD 50 ' 7 collabora- 
tions. Their results are consistent with each 
other and slightly lower than the two-flavor 
data. My average is m s (2 GeV) = 78 ± 
10 MeV. 

Determination of light quark mass m or 
the ratio m s /fh is sensitive to the chiral ex- 
trapolation. At the leading order of m q the 
ratio is given by the physical meson masses as 
to^/to 2 - 1 = 25.9, and a NLO ChPT analy- 
sis yields 24.4±1.5 51 . The lattice calculation 
can be used to improve this estimate. The 
small quark mass reached by the MILC sim- 
ulation enabled them to include NLO ChPT 
terms in the fit 50 ' 7 as well as the correc- 
tion terms to describe the taste symmetry 
breaking 52 and higher order effects. Their 
result 27.4±4.2 is consistent with the NLO 
ChPT analysis but slightly higher, suggest- 
ing non-negligible higher order effect. 

3.3 Heavy quark masses 

The charm quark is not too heavy to describe 
with the 0(a)-improved Wilson fermion ac- 
tion adopting the naive estimate of discretiza- 
tion effect 0((am c ) 2 ). It can in princi- 
ple be eliminated by taking the continuum 
limit, which is feasible in the quenched ap- 
proximation and precise results to c (to c ) = 
1.30(3) 53 and 1.32(3) 54 GeV are obtained. 
If one takes the non-relativistic dynamics 
of charm quark inside the D( s ) meson into 
account 55,56 , the discretization error is not 
as large as 0((am c ) 2 ). Recent work indi- 
cates that the discretization effect is much 
smaller , and an unquenched calculation is 



d In these works, the matching is done perturbatively. 
The QCDSF-UKQCD collaboration 47 calculated the 
non-perturbative matching factor for the VWI de- 
termination, and found it larger by about 20% than 
the one-loop estimate. The central values of 36,6 are 
taken from the AWI definition, however. 



being performed. 

For the bottom quark the conventional 
approach fails for the lattice scale 1/a ~ 
2-3 GeV. Instead, the heavy quark effec- 
tive theory (HQET) is a good approximation 
up to corrections of order 0(Aq CD /to;,) ~ 
30 MeV. Higher order perturbation theory is 
essential for the matching of nib in order to 
avoid large corrections due to power diver- 
gences. The two-loop calculation was done 
sometime ago 58 and the three-loop calcula- 
tion has been performed recently 59 , reduc- 
ing the error to the 40 MeV level. Avail- 
able two-flavor QCD calculations combined 
with the two-loop matching yield fhb{fhb) = 
4.21(7) MeV 60 ' 59 and 4.25(11) MeV 61 . For 
the latter, carefully estimated uncertainties 
in the lattice scale and strange quark mass 
dominate the error bar, which is expected to 
be reduced by 2+1-flavor calculations. 

Recently, a non-perturbative method to 
match HQET onto QCD has been formulated 
and tested on quenched lattices 62 . Another 
method to calculate b quark mass without re- 
course to HQET has also been proposed 54 . 
These methods may enable us to further re- 
duce the systematic error. 

4 Kaon physics 

4-1 Determination of \V US \ 

The best known method to determine \V US \, 
or the Cabibbo angle, is to use the semi- 
leptonic Ki% decays. The relevant form fac- 
tor /+(0) is normalized to one in the SU(3) 
limit (to = m s ) , and the correction starts at 
the second order in m s — m 63 . Calculation 
of the correction in a quark model yielded 
/+(0) = 0.961 (8) 64 . Further improvement re- 
quires non-perturbative method to calculate 
/+(0), and first quenched lattice calculation 
has been done recently 65 using double ratios 
as in the \ V c b\ calculation 66 ' 67 . They reported 
0.960(5)(7). 

V us | can also be determined through the 
leptonic decay AT ± — > (ir^Vu, once the decay 
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Figure 3. Quenched lattice results for Bk (2 GeV) 
versus lattice spacing a. Results with stag- 
gered fermion (JLQCD 69 , Lee et al. 70 , Gamiz 
et al. 71 ), Wilson- type fermions (SPQcdR clover 72 , 
SPQcdR Wilson 73 ) twisted mass fermion 74 , domain- 
wall fermion, (CP-PACS 45 , RBC 75 ' 76 ), and overlap 
fermion (DeGrand 43 , Garron et al. 77 ). Horizontal 
lines show my average. 



constant Jk is known theoretically. This has 
been attempted 68 using the recent MILC re- 
sult f K /fn = 1. 210(4) (13) 7 . It is notable 
that the error is now comparable to the semi- 
leptonic determination, i.e. ~ 1%, and the 
result for |V^ S | is consistent. 



4-2 Kaon B parameter 

Calculation of Bk is one of the major goals 
of lattice QCD calculations. Bk is a matrix 
element of a AS = 2 four-quark operator of 
VV + AA chiral structure sandwiched by K 
and anti-if states. In the lattice calculation 
its chiral structure must be maintained in or- 
der to avoid large contamination from wrong 
chirality operators. Therefore, most lattice 
calculations have been done using fermion 
formulations which respect chiral symmetry. 

In the quenched approximation, the 
"benchmark" result was given by the JLQCD 
collaboration 69 in 1997 using (unimproved) 
staggered fermion: B^ /=0 (2GeV) = 0.63(4). 
As shown in Figure 3 they calculated Bk 
at several values of lattice spacing down to 



~ 0.04 fm and extrapolated to the contin- 
uum limit assuming theoretically expected 
discretization effects. More recent improved 
staggered fermion results 70 ' 71 are consistent 
with the JLQCD's continuum limit already 
at large lattice spacings a ~ 0.1-0.2 fm. 

In the past few years, many calculations 
have been performed using the Ginsparg- 
Wilson fermions (domain-wall and overlap) 
45,75,76,43,77_ tj,^ resu ]^ s are shown by tri- 
angles in Figure 3. They are slightly lower 
than the JLQCD's continuum limit, and dis- 
cretization error is substantially reduced. For 
some of these, the renormalization factor is 
computed non-perturbatively. 

Another method to protect the lattice 
operator from operator mixing is provided by 
chirally twisted mass lattice QCD 78 . The nu- 
merical result for Bk 74 is almost flat in a and 
consistent with the Ginsparg- Wilson fermion 
results 6 . 

In view of these improved calculations, 
which have smaller discretization effect and 
non-perturbative renormalization, I recom- 
mend a slightly lower value for a quenched 
world average, B^ /= °(2GeV) = 0.58(4), 
which is shown in Figure 3 by horizontal lines. 

Dynamical quark effects were previ- 
ously estimated from unimproved stag- 
gered simulations 80 ' 81 ' 82 ' 83 as only slightly 
negative 84 or positive 85 . Due to the large 
discretization errors inherent in the unim- 
proved staggered fermion, it was difficult to 
disentangle the sea quark effect from the dis- 
cretization effect. This year, new unquenched 
calculations have appeared using the 0(a)- 
improved Wilson fermion 86 , improved stag- 
gered fermion 71 , and domain- wall fermion 87 . 

Figure 4 shows sea quark mass (or mass 
of pion composed by sea quarks) dependence 
of Bk- The RBC calculation with dynami- 
cal domain-wall fermion 87 is very precise and 



e A related proposal has been made for Wilson 
fermions 79 to avoid the wrong chirality operators. A 
numeical result 73 suggests large O(a) discretization 
effect for unimporved Wilson fermion. 
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Figure 4. Unquenched B K . Data from UKQCD sb 
and RBC 87 are shown as a function of pion mass 
squared. The right panel shows quenched world av- 
erage. 



shows a tread to decrease toward the chiral 
limit. The SU(3) breaking m s ^ ma effect is 
estimated to be about —3%. Since the result 
is still preliminary and the slope in the chiral 
extrapolation seems to rely on a single point 
at rn 2 PS ~ 0.25 GeV 2 , I take this as an indi- 
cation of quenching error rather than a cen- 
tral value of unquenched QCD, and recom- 
mend an average B K (2GeV) = 0.58(4)(±°), 
where the second asymmetric error reflects 
the quenching error. 



5 Heavy quarks 

5.1 D meson decays 

The CLEO-c and BES III experiments 
promise to measure the D( s ) decays at a few 
% accuracy, which provides a stringent test of 
lattice simulation of heavy quarks. The most 
relevant hadronic quantities to be calculated 
on the lattice are the leptonic decay constants 
and semi-leptonic decay form factors. 

As in the charm quark mass calculation 
the large discretization effect of 0((am c ) 2 ) 
must be eliminated for precise calculation of 
meson decay constants, unless one uses 
effective theory approaches. Continuum ex- 
trapolation has recently been performed in 
the quenched approximation, yielding fu s = 
252(9) MeV 88 and 240(5)(5) MeV 89 . They 



are consistent with the previous world aver- 
age 230(14) MeV 90 from the calculations us- 
ing the Fermilab heavy quark formulation 55 , 
which applies the idea of HQET for the 
Wilson-type lattice fermion. 

The Fermilab-MILC-HPQCD collabora- 
tion has carried out calculations of fu and 
fu. in the presence of 2+1-flavors of stag- 
gered sea quarks 91 . Their result is fr>, = 
263(tg)(24) MeV. The first error is statis- 
tical, and the second error reflects their es- 
timate of systematic error due to matching 
of the heavy quark action and current to 
QCD (7%), discretization effects of the light 
quark (4%), charm quark mass determina- 
tion (4%), etc. To improve the heavy quark 
matching, one requires perturbative calcula- 
tion for complicated lattice actions including 
many 0{l/m,Q) terms, and such work is in 
progress using automated perturbative cal- 
culation technique 92 . Other errors can be 
reduced by increasing the computing power, 
and it would be feasible to reduce the total 
error to the 5% level. 

The semi-leptonic decay D — > ir£v, Klv 
form factors have also been calculated by the 
same group, and the results are fS^^ = 
0.64(3)(6) and ff^ K = 0.73(3)(7), which 
can be used to determine |V CS | and \V c d\ 
with experimental inputs from BES 93 and 
CLEO 94 . Such determination is currently 
consistent with the CKM unitarity and the 
error is dominated by lattice calculation. If 
we can reduce the error to the 5% level or 
better, it will provide interesting test of the 
CKM unitarity. 



5.2 B meson mixings 

The mass difference in the B° — B° system 
A Aid is proportional to f B BB\V t d\ 2 , and thus 
gives a constraint on the CKM element \Vtd\, 
provided that the corresponding hadronic 
matrix elements, B meson decay constant fs 
and B parameter Bb, are theoretically calcu- 
lated. The analogous mass difference in the 



For Publisher's use 



rt=b 

AliKhan et al. (1998) 
JLQCD (1999) 
CPTACS (2001) 



JLQCD (1998) I — e— 

EI-Khadra el al. 11998) j-e 1 

MILC Q 998) h- 8 1 



CP-PACS (2000) 
MILC (2002) 

Becirevic et al. (1998) 
UKQCD (2000) 
Lellouch-Lin (2000) 
Becirevic et at (2000) 

de Divitiis et al. (2003) 
ALPHA (2003) 

N J= 2 

Collins etal. (1999) 
CP-PACS (2000) 
CP-PACS (2001) 
MILC (2002) 
JLQCD (2003) 

N = 2+1 



ingate et al. (2003) 
I ... I 



Wii 



120 160 200 240 
f d [MeV] 



280 



Figure 5. Lattice calculations of /s s - Quenched (7V^ 
= 0) results are obtained with the NRQCD action 
(Ali Khan et al. 95 , JLQCD 96 , CP-PACS 97 ), Fer- 
milab formulation (JLQCD 98 , El-Khadra et al. 99 , 
MILC 100 , CP-PACS 101 , MILC 102 ), and conventional 
Wilson fermions (Becirevic et al. 103 , UKQCD 104 , 
Lellouch-Lin 105 , Becirevic et al. 106 , de Divitiis et 
al. 89 , ALPHA 107 ). Unquenched (N f = 2, 2+1) cal- 
culations are Collins 108 , CP-PACS 101 ' 97 , MILC 102 , 
JLQCD 109 , Wingate et al. 110 . 



— B® mixing AM S can be used to reduce 
the theoretical uncertainty by considering a 



ratio AM s /AM d = e(M B JM B )\V ts /V t 



where a ratio £ = (fB a Bg 2 )/(fBB B /z ) to de- 
scribe the SU(3) breaking effect is the quan- 
tity to be estimated theoretically. 

In the lattice calculation, the difficulty 
to treat heavy quarks on the lattice has been 
essentially solved by introducing the HQET 
based lattice actions, and the results with dif- 
ferent formulations are in good agreement in 
the quenched approximation within the sys- 
tematic uncertainty of order of 15%, as shown 
in Figure 5 for fg s . 

Unquenched simulations have been done 
since 1999108,101,97,102 p igure 5 shows the 

results for fg s , which is relatively insensi- 
tive to the problem of chiral log and the chi- 
ral extrapolation is more reliable. The un- 
quenched results seemed slightly higher than 
the quenched results, but it depends on which 
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Figure 6. A ratio (/b s -M^/ 2 )/(/b Mg^ 2 ) as a func- 
tion of light and strange quark mass ratio. Lattice 
data are from JLQCD 109 and HPQCD 116 . 



quantity is used to set the lattice scale, e.g. 
if we set the scale with the heavy quark po- 
tential (rrj) the quenched and unquenched re- 
sults are consistent with each other. On the 
other hand, the disagreement between the 
most recent two calculations, JLQCD (Nt = 
2) 109 and Wingate et al. (N f = 2 + l) 110 , re- 
mains even if one uses the common scale such 
as ro. At present, it is not clear if it comes 
from the dynamical strange quark effect or 
from other (underestimated) systematic er- 
rors. For an average of unquenched calcula- 
tions I recommend Jb s = 230 ± 30 MeV. 

For the heavy-light meson decay con- 
stant, ChPT predicts a non-analytic quark 
mass dependence 111 oc (1 + 3g 2 )m 2 lnm 2 , 
whose effect could become important in the 
chiral extrapolation of fs ■ The coupling g de- 
scribe the HH*ir interaction, such as BB*ir 
or DD*tt. For the DD*tt interaction, there is 
an experimental measurement g — 0.59(7) 112 , 
and quenched lattice calculations 113,114 are 
also consistent with this. 

In Figure 6 the chiral extrapolation is 
shown for a ratio (/q, )/(/b^b ) as 
a function of m q /m s f. It goes up to- 
ward the physical point m q /m s « 1/25. 
The results from the JLQCD collaboration 109 



above m q /m s 



0.5 are consistent with 



■''The plot is an update of that shown by Kronfeld at 
Lattice 2003 115 . 
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Table 1. My averages of B mixing parameters com- 
pared to those by Lellouch at ICHEP 2002 2 . f B 



and f E 



B 1/2 



are given in units of MeV. 



a linear dependence, but the chiral log 
could raise the chiral limit as shown by the 
solid curve. Recently the data from the 
HPQCD collaboration 116 with the staggered 
sea quarks become available at smaller quark 
masses, which are consistent with the ex- 
pected chiral log behavior, although more 
statistics is needed to be conclusive. I fit the 
both data with a function including the chiral 
log and obtain f B Jf B = l-22(+|). 

The large uncertainty due to the chiral 
log can be avoided by constructing a dou- 
ble ratio 117 (f B J f B )/{f D J fn), in which the 
chiral log term cancels at the leading order 
of 1/M. The JLQCD collaboration calcu- 
lated this ratio and obtained a very precise 
result 1.01 (l) 118 . Once the D meson decay 
constants are measured precisely at CLEO-c 
and BES III, this strongly constrains the ra- 
tio f Bs /f B . One could also consider another 
double ratio {f B Jf B )/{fK/U), for which 
the effect of chiral log is numerically small 119 . 

Fortunately, the coefficient of the chi- 
ral log term is small for the B parame- 
ter B B , and the unquenched lattice calcu- 
lation shows no significant sea quark mass 
dependence. Therefore, the results of the 
JLQCD collaboration 109 will be robust as far 
as the chiral extrapolation is concerned. Un- 
quenched results are B B {rrib) — 0.84(3) (6) 
and B B JB B = 1.02(2) (t 6 ), which are con- 
sistent with the previous quenched results 120 . 



UKQCD (1999) 
Abada et al. (2000) 
> El-Khadra et al. (2001) 
< JLQCD (2001) 



N=2 + l 




Fermilab (2004) 
HPQCD (2004) 



Figure 7. The B — + nlv decay form factors fj r (q 2 ) 
(upper cluster of points) and foil 2 ) (lower clus- 
ter) from quenched and unquenched lattice QCD. 
Quenched (N f = 0) results are from UKQCD 123 , 
Abada et al. 124 , El-Khadra et al. 125 , JLQCD 126 , 
while the unquenched results are from Fermilab 121 
and HPQCD 122 . 



My averages for the B — B mixing param- 
eters are summarized in Table 1. The aver- 
ages at ICHEP 2002 provided by Lellouch 2 
are also listed for comparison. The main dif- 
ference is the treatment of the uncertainty 
due to the chiral extrapolation. Because an 
indication of the chiral log is already found 
I take the central value from the fit includ- 
ing the chiral log, while it was given in the 
asymmetric second error in 2 . 

5.3 B meson decays 

The B meson semi-leptonic decay form fac- 
tors are needed in the determination of \V u b\ 
and \Vcb\ from exclusive decay modes B — » 
■kIv and B -> D { *Hv. 

The first unquenched calculations of the 
B — > ~kIv form factors have recently been 
presented by the Fermilab lattice 121 and 
the HPQCD 122 collaborations, which both 
worked on 2+l-flavor gauge field ensembles 
produced by the MILC collaboration. Fig- 
ure 7 shows their results together with the 
previous quenched calculations 123 ' 124 ' 125 ' 126 . 
We do not observe any significant difference 
between quenched and unquenched calcula- 
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tions within relatively large statistical errors. 

Since lattice calculations are possible 
only in the high q 2 region where the recoil 
momentum of the daughter pion is small, the 
extraction of \V u b\ using these lattice results 
requires experimental data in the same kine- 
matical region. Such analysis was done by 
the CLEO collaboration 127 . They provided 
the data of the decay rate in three q 2 bins, 
and the highest bin (above 16 GeV 2 ) cor- 
responds to the region of the lattice calcu- 
lation. The result using an average of the 
four quenched lattice calculations \V u b\ = 
(2.88 ± 0.55 ± 0.30 ±±IH ±0.18) x 10~ 3 still 
has large errors (due to statistical, system- 
atic, lattice, and piv form factor dependence, 
in the order given), but we expect better mea- 
surements from BaBar and Belle in near fu- 
ture. In fact, a preliminary result with higher 
purity (and thus with smaller experimental 
systematic error) was reported at this con- 
ference by the Belle collaboration 128 . 

The Fermilab lattice collaboration has 
updated 121 the calculation of B — > Dtv form 
factor at zero recoil, which can be used in 
the precise determination of | V c b \ ■ Their new 
result with 2+1 flavors of dynamical quarks, 
J~b^d{^) = 1.075(18)(15), is consistent with 
the previous quenched calculation 66 . Similar 
calculation of the B — > D^*'£u form factor 67 
is in progress. 

6 Conclusions 

The most important progress in lattice QCD 
calculations in the past few years is the in- 
clusion of sea quark effects. By performing 
the real simulations with dynamical quarks, 
the unknown errors due to quenching in the 
previous calculations are being eliminated for 
many of the quantities discussed here. 

In unqucnchcd calculations, the non- 
analytic quark mass dependence appearing 
from the pion loops must be taken into ac- 
count in the chiral extrapolation. Without 
enough data points below mps — 500 MeV, 




Figure 8. Constraints on the CKM unitarity triangle 

1 /2 

using the most recent lattice results for Bk , JB a BJ , 
and £. Plot is from the UTfit collaboration 129 . 



the extrapolation induces large systematic 
uncertainty, since we do not know where the 
chiral log effect becomes important. In this 
respect the advantage of staggered fermion 
is clear: the simulation is the fastest and 
therefore one can reach the chiral regime. 
Employing the improved staggered fermions, 
the HPQCD-MILC-UKQCD-Fermilab group 
have recently presented several results with 
2+1 flavors of dynamical quarks. 

The drawback of the staggered fermion 
is its complicated taste structure, in particu- 
lar the fourth-root trick introduced to repre- 
sent single dynamical flavor. It is not proved 
that the fourth-root of the staggered fermion 
determinant is the same as a local quantum 
field theory. Such a proof is essential, as it 
would provide a theoretical basis that these 
calculations correspond to the real QCD. Un- 
til such a proof is available, the simulation 
with other fermion formulations without such 
taste structure should be pursued. 

The constraints on the CKM unitarity 
triangle using the recent lattice results are 
shown in Figure 8. Although the error bars 
are not significantly reduced over the last sev- 
eral years, the uncertainty due to quenching 
is now basically eliminated, and the calcula- 
tions are from the first principles. Much more 
work is needed to reduce other systematic er- 
rors, and the ideas are already being tested 
within the quenched approximation. 
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